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1.0  Problem  Statement 

This  report  summarizes  the  work  performed  under  the  contract  DAAL03-91-C-0052, 
“Modeling  of  Chaotic  Ructuating  Channels  for  Signal  Detection.”  This  work  was  per¬ 
formed  from  July  1991  to  April  1995. 

The  detection  of  signals  propagated  through  a  fluctuating  medium  is  a  significantly  chal¬ 
lenging  problem  and  is  of  significant  practical  interest.  The  objective  of  this  research  was 
to  investigate  the  tools  and  techniques  used  to  model  fluctuating  propagation  as  nonlinear 
systems  and  to  analyze  propagation  data  in  the  context  of  a  nonlinear  systems  model. 

The  objectives  of  this  program  were  as  follows: 

•  Develop  signal  processing  algorithms  and  software  applicable  to  the  modeling  and 
analysis  of  nonlinear  system.  We  were  to  develop  tools  and  techniques  which  had  appli¬ 
cability  to  a  range  of  nonlinear  processing  problems  and  to  apply  them  both  to  the  prob¬ 
lem  at  hand,  propagation  of  signals  through  a  fluctuating  medium,  as  well  as  to  other 
problem  which  seemed  to  be  amenable  to  these  techniques. 

•  Develop  signal  detection  algorithms  which  are  well  suited  to  the  exploitation  of  nonlin¬ 
ear  propagation  and  interference  phenomenon.  We  were  to  develop  signal  processing 
algorithms  which  would  allow  us  to  exploit  the  nonlinear  description  of  systems  to 
extract  information  about  signals  of  interest. 

•  Perform  analyses  of  acoustic  signals  propagated  through  the  atmosphere  to  determine  if 
the  propagation  effects  can  be  well  modeled  as  a  nonlinear  system.  We  were  to  design  a 
propagation  experiment  to  collect  air  acoustics  data  to  determine  if  propagation  of  sig¬ 
nals  in  the  real  world  could  be  modeled  as  the  response  of  a  nonlinear  system. 

•  Identify  new  opportunities  for  the  application  of  nonlinear  processing  to  additional 
areas  of  interest.  We  were  to  investigate  datasets  of  interest  to  identify  new  opportuni¬ 
ties  for  the  application  of  nonlinear  modeling  and  processing. 

2.0  Summary  of  Significant  Results 

During  the  course  of  this  project  we  achieved  significant  progress  against  the  four  main 
objectives.  Most  of  this  work  has  been  documented  in  published  material  and  will  be  only 
summarized  here. 

2.1  Algorithms  and  Software  for  the  Modeling  and  Analysis  of 
Nonlinear  Systems 

As  part  of  our  contract  work  we  have  developed  a  wide  range  of  algorithms  for  the  analy¬ 
sis  and  modeling  of  nonlinear  systems  [3, 4, 8, 11, 12, 14, 22, 24, 41].  These  tools  include 
algorithms  for  the  following: 

•  establishing  the  proper  method  data  selection  for  modeling  a  data  stream  as  the  output 
of  a  nonlinear  system; 
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•  determining  the  proper  embedding  dimension  of  data  for  analysis  as  the  output  of  a 
nonlinear  system; 

.  quantitatively  differentiating  between  the  output  of  low  dimensional  nonlinear  system 
and  “noise”; 

.  determining  invariant  characteristics,  e.g.  Lyapunov  exponents,  of  a  nonlinear  system 
from  its  outputs;  and 

•  building  a  predictive  model  of  a  nonlinear  system  which  can  be  used  for  data  extrapola¬ 
tion  or  for  control  of  a  nonlinear  system. 

The  software  for  these  algorithms  has  been  made  available  to  researchers  worldwide 
through  the  University  of  California,  San  Diego  and  is  being  incorporated  into  “tool¬ 
boxes”  for  use  in  standard  signal  processing  packages  such  as  Matlab. 


2.2  Signal  Detection  Algorithms 

As  part  of  the  contract  work  we  developed  algorithms  for  the  separation  of  signals  pro¬ 
duced  by  low  dimensional  nonlinear  systems  from  signal  produced  by  noise  sources  [1,2, 
22, 24].  This  work  resulting  in  the  development  of  algorithms  which  successfully  merged 
clisical  signal  processing  techniques,  such  as  the  Cramer-Rao  lower  bound  and  the  Vit- 
erbi  algorithm,  with  nonlinear  systems.  With  this  work  we  developed  the  following  signif¬ 
icant  results: 

•  Demonstration  that  the  noise  reduction  potential  for  nonlinear  systems  is  globally  the 
same  as  the  noise  reduction  potential  for  linear  systems  but  the  local  properties  of  noise 
reduction  in  nonlinear  systems  is  drastically  different. 

•  The  signal  structure  of  a  chaotic  signal  can  be  exploited  by  a  signal  detection  algorithm 
which  has  been  trained  to  the  chaotic  system  of  interest  to  provide  a  significant  process¬ 
ing  gain  compared  to  a  conventional  colored  noise  detection  system. 


2.3  Analysis  of  the  Propagation  of  Acoustic  Signals 

As  part  of  our  contract  work  we  collected  and  analyzed  acoustic  propagation  data  from  a 
real-world  environment.  The  data  collection  experiment  took  place  at  the  US  Army  White 
Sands  Missile  Range  in  White  Sands,  New  Mexico  in  December  1992.  In  preparation  for 
this  experiment  we  defined  a  wide  range  of  signals  of  interest  for  propagation  through  the 
atmosphere;  we  determined  the  proper  sensor  locations;  we  defined  the  data  collection 
hardware  and  software  that  would  be  used;  we  coordinated  with  the  Atmospherics  Sci¬ 
ences  Laboratory  at  White  Sands  to  use  some  of  their  equipment;  and  we  coordinated  with 
another  Sanders  contract,  this  one  to  ARDEC,  to  jointly  collect  data  (Appendix  A). 

The  data  was  collected  at  White  Sands  over  a  period  of  three  day  in  the  morning,  at  mid¬ 
day,  and  in  the  evening  for  a  total  of  ten  separate  data  collections.  Each  collection  con¬ 
tained  sixteen  channels  of  recorded  data  with  each  recording  lasting  between  ten  and  thirty 
minutes.  The  sensor  locations  included  an  array  of  microphones  on  the  ground  as  well  as 
microphones  placed  on  a  tower  at  various  elevations  off  the  ground. 
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The  collected  was  analyzed  through  a  number  of  nonlinear  processing  techniques  to 
determine  if  any  significant  nonlinear  structure  could  be  found  in  the  data  (Appendix  B). 
We  were  not  able  to  identify  any  significant  low  dimensional  nonlinear  behavior  in  the 
Hata  Additionally,  classical  signal  processing  techniques  were  used  to  analyze  the  data 
and  these  demonstrated  that  the  data  was  best  modeled  as  simple  direct  path  channel  with 
random  fluctuations  on  the  propagation.  The  data  is  dominated  by  acoustic  interfere  due  to 
acoustic  sources  in  the  area  over  which  we  had  no  control  and  by  wind  noise  which  signif¬ 
icantly  interferes  with  the  measurement  of  the  propagated  signals. 


2.4  Identification  of  Other  Applications  for  Nonlinear  Processing 

Because  of  the  lack  of  significant  nonlinear  behavior  in  the  White  Sands  data,  we 
expanded  our  data  analysis  work  to  include  additional  datasets.  In  this  portion  of  the  work 

we: 

•  examined  wind  tunnel  data  collected  on  a  jeep  at  the  NASA  Langley  wind  tunnel, 

•  examined  data  collected  by  the  Navy  on  a  buoyant  test  vehicle  [19], 

•  analyzed  various  climatic,  laser,  and  fluidized  bed  data  [various];  and 

•  developed  an  experimental  plan  to  provide  high  quality  air  acoustic  flow  noise  data 
[46]. 

The  results  of  the  wind  tunnel  data  indicated  that,  at  least  for  some  sensors  on  the  vehicle 
at  some  wind  speeds,  that  nonlinear  behavior  was  observed.  Unfortunately  the  data  ana¬ 
lyzed  here  was  collected  under  conditions  which  were  not  completely  defined  by  the  data 
we  received.  The  buoyant  vehicle  experiment  was  much  better  instrumented  and  con¬ 
trolled.  The  results  of  the  buoyant  test  vehicle  experiment  data  clearly  showed  the  transi¬ 
tion  from  laminar  to  turbulent  flow  and  regions  of  low  dimensional  behavior  were 
observed.  We  took  the  work  done  here  and  extrapolated  to  define  an  experiment  that  could 
be  performed  for  the  propagation  of  air  over  a  surface  in  a  controlled  environment  and 
with  a  high  resolution  data  acquisition  system  to  develop  a  high  quality  data  set  which  will 
demonstrate  the  transition  from  laminar  to  turbulent  flow  in  air  propagation  and  which 
may  lead  to  insight  into  the  mechanisms  by  which  acoustical  noise  is  generated  by  the 
wind. 
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Chaotic  Laser  Dynamics”,  submitted  to  Phys.  Rev.  Letters,  July,  1994. 

[39]  Abarbanel,  H.  D.  I.,  Z.  Gills,  R.  Roy,  and  C.  Liu,  “Quantum  Fluctuations  in  the  Oper¬ 
ation  of  a  Chaotic  Laser”,  submitted  to  Physical  Review  E,  January,  1995. 

[40]  Abarbanel,  H.  D.  L,  P.  F.  Rowat,  and  A.  I.  Selverston,  “Nonlinear  Dynamics  of  the 
Lobster  Gastric  Mill  Central  Pattern  Generator”,  to  be  submitted  Fall,  1994. 

[41]  Kennel,  M.  B.  and  H.  D.  I.  Abarbanel,  “False  Strands  and  False  Neighbors”,  to  be 
submitted  to  Physical  Review  E,  Winter,  1995. 

[42]  Abarbanel,  H.  D.  L,  N.  F.  Rul'kov,  and  M.  M.  Sushchik,  “Blending  Chaotic  Attractors 
Using  Chaotic  Synchronization”,  submitted  to  Physical  Review  Letters,  October,  1994. 
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[43]  Abarbanel,  H.  D.  I.,  T.  Prison,  and  L.  Sh.  Tsimring,  “Chaotic  Methods  for  Image  Pro¬ 
cessing”,  ICASSP,  Detroit,  MI,  May,  1995. 

[44]  Abarbanel,  H.  D.  I.,  T.  Prison,  R.  Katz,  and  J.  Cembrola,  “Nonlinear  Analysis  of 
Environmental  Distortions  of  Continuous  Wave  Signals  in  the  Ocean  ,  submitted  to  J. 
Acoust.  Soc-  Am,,  November,  1994. 

[45]  Abarbanel,  H.  D.  I.  and  Upmanu  Lall,  “Nonlinear  Dynamics  of  the  Great  Salt  Lake: 
System  Identification  and  Prediction”,  submitted  to  J.  of  Climate,  November,  1994. 

The  following  material  represents  unpublished  reports  developed  under  this  contract: 

[46]  Dingus,  P,  Lang,  S.,  Myers,  C.,  Nation,  R.,  Abarbanel,  H.,  “Basic  Physics  and  Appli¬ 
cations  of  Nonlinear  Dynamics,”  submitted  to  ARO,  April  1995. 

[47]  Abarbanel,  H.  D.  I.,  “Eikonal  approximation  techniques  in  elastic  scattering  and  pro¬ 
duction  processes”.  Lecture  Notes  in  Physics:  Strong  Interaction  Physics,  (Summer 
School  Lectures  at  the  Canadian  School  of  Physics)  (1973). 

[48]  Abarbanel,  H.  D.  I.,  Gribov,  V.N.,  and  Kanchelli,  O.V.,  “Couplings  of  the  vacuum  tra- 
jectory.“ 

[49]  Abarbanel,  H.  D.  I.,  “Strongly  mixed  stationary  turbulence”,  LBL  Report  (1979). 

[50]  Abarbanel,  H.  D.  I.,  Grebogi  C.,  and  Kaufman,  A.,  “Statistical  analysis  of  a  determin¬ 
istic  orbit”,  Proc.  Inti.  Conf  on  Plasma  Physics  (1980). 

[51]  Abarbanel,  H.  D.  I.,  Grebogi,  C.,  and  Kaufman,  A.  “Statistical  analysis  of  a  determin¬ 
istic  mapping”.  Report  to  the  Plasma  Physics  Division  of  the  APS,  (1980). 

[52]  Abarbanel,  H.  D.  I.,  Grebogi,  C.,  and  Kaufman,  A.,  “Resonance  diffusion  in  an 
eikonal  field”.  Report  to  the  Plasma  Physics  Division  of  the  APS,  (1980). 

[53]  Abarbanel,  H.  D.  I.,  “Irreversibility  and  transport  equations  in  nonlinear  deterministic 
dynamics”.  Seventh  West  Coast  Statistical  Mechanics  Conference,  LBL-12683,  (1981). 

[54]  Abarbanel,  H.  D.  I.,  and  Crawford,  J.D.,  “Bifurcation  theory  and  plasma  instabili¬ 
ties”,  Report  to  the  Plasma  Physics  Division  of  the  American  Physical  Society  (1982). 

[55]  Abarbanel,  H.  D.  1.,  Brown,  R.,  and  Yang,  Y.M.,  “Equations  of  State  with  Maximum 
Symmetry  in  Inviscid  Fluids”,  (October,  1987). 

[56]  Abarbanel,  H.  D.  I.,  Brown,  R.,  and  Yang,  Y.M.,  “Symmetries  and  Conserved  Quanti¬ 
ties  for  Inviscid,  Compressible  Fluids  in  Two  and  Three  Space  Dimensions”. 

[57]  Abarbanel,  H.  D.  I.,  Brown,  R.,  and  Yang,  Y.M.,  “Symmetries  and  Conserved  Quanti¬ 
ties  for  Inviscid,  Compressible  Fluids  in  One  Space  Dimension”. 
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[58]  Abarbanel,  H.  D.  1.  and  Hjorth,  R,  “Surface  Instability  of  a  Driven  Non-Newtonian 
Film  Flow”. 

[59]  Abarbanel,  H.  D.  I.  and  J.  B.  Kadtke,  “Information  Theoretic  Methods  for  Determin¬ 
ing  the  Minimum  Embedding  Dimension  for  Strange  Attractors  . 

The  following  represent  material  which  is  still  under  preparation  but  which  was  begun 
under  this  contract: 

[60]  Abarbanel,  H.  D.  1.,  Analysis  of  Observed  Chaotic  Data,  to  be  published  by  Springer- 
Verlag,  Summer,  1995. 

[61]  Abarbanel,  H.  D.  L,  M.  I.  Rabinovich,  and  Lev  Sh.  Tsimring,  Spatio-temporal  Chaos: 
Structures  and  Experimental  Analysis  to  be  published  by  Springer- Verlag. 

[62]  Richard,  M.  D.  and  H.  D.  I.  Abarbanel,  “Estimating  in  Chaos:  Cramer-Rao  Bounds 
for  Chaotic  Tims  Series”,  UCSD/MIT  Preprint,  Spring,  1993.  To  be  submitted  to  Journal 
of  Nonlinear  Science . 


4.0  List  of  Participating  Scientific  Personnel 

The  work  reported  on  in  this  contract  was  performed  in  cooperation  of  Sanders,  under  the 
direction  of  Dr.  Cory  Myers  and  at  the  University  of  California,  San  Diego  under  the 
direction  of  Professor  Henry  Abarbanel. 

The  following  personnel  from  Sanders  participated  on  this  project: 

Dr.  Cory  Myers 
Dr.  Steven  Lang 
Mr.  Edward  Real 
Ms.  Frances  Shin 

Mr.  Andrew  Singer  (summer  smdent  from  MIT) 

Mr.  Jeffery  Ludwig  (summer  student  from  MIT) 

The  following  personnel  from  the  University  of  California,  San  Diego  participated  on  this 
project: 

Professor  Henry  Abarbanel 

Dr.  Matt  Kennel  (received  his  Ph.D.  during  the  course  of  this  project) 

Dr.  J.  J.  (“Sid”)  Sidorowich 
Dr.  L.  Sh.  Tsimring 
Dr.  Reggie  Brown 

Dr.  Mischa  Sushchik  (will  receive  his  Ph.D.  this  fall) 


5.0  Inventions 

No  inventions  were  patented  under  this  contract. 
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6.0  Bibliography 

See  the  list  of  publications  and  technical  reports  in  section  2. 


7.0  Appendices 

Enclosed  are  copies  of  material  which  was  developed  under  this  contract  but  which  has 
not  been  previously  published  or  submitted. 

Appendix  A  -  Experimental  Plans  -  Documents  describing  the  experimental  setup  for  the 
White  Sands  data  collection  of  December  1992.  Included  here  is  the  associated  data  col¬ 
lection  plan  for  the  ARDEC  data  collection  which  shared  facilities  and  equipment  with  our 
data  collection. 

Appendix  B  -  Sanders  Internal  Data  Analysis  Reports  -  A  collection  of  data  analysis 
reports  relating  to  the  analysis  of  the  White  Sands  data  and  to  the  NASA  wind  tunnel  data. 
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Appendix  A 


ARO  Test  Plan 


Lockheed  Sanders,  Inc. 


Modeling  of  Chaotic  Fluctuating  Channels  for  Signal 

Detection 

Experiment  Description 


1  Objective 

The  obiective  of  the  proposed  data  collection  is  to  provide  high  quality  data  sets  that  can  be  ana- 

Izedlr^work  for  aL  under  contractDAAL03-91-^^^^^ 

to  apply  recent  developments  in  the  area  of  nonlinear  systems  to  the  problem  of  detecting  acoustic 
sirrnds^that  have  been  propagated  through  the  atmosphere.  Our  objectives  are  to  extend  our 
Sder^lnl '  of  nonliL^  systems,  to  apply  the  theory  of  nonlinear  systems  to  the  modeling  of 
fluctuations  in  the  atmosphere  and  their  effects  on  the  propagation  of  signals,  and  to  develop  pro 
cessing  approaches  for  the  improved  detection  of  signals  as  propagated  through  fluctuating 

nels. 

The  program  is  been  performed  jointly  by  .he  Signal  Processing  Center  TKh“do^ 

(SPCOT)  ar  Lockheed  Sanders,  Inc.  and  die  Insdrure  for  Nonlinear  Scrence  ONLS)  “ 
sitv  of  California,  San  Diego  (UCSD).  SPCOT  has  extensive  expenence  in  the  processing  o 
acoustic  signals,  panicularly  for  detecnon  and  classification.  INLS  has  extensive  expenence  in 
the  theory  of  nonlinear  systems  and  in  applying  these  models  to  physical  systems. 

The  particular  objectives  of  the  experiment  are  to  provide  high  quality  acoustic  propagation  data 
serwelsire  dat  that  can  be  used  to  characterize  die  fluctuations  in  die  atmosphere,  how  these 
flucmldons  propagate,  and  how  the  dansmission  of  signals  are  affected  by  these  flucmatrons. 

2  Approach 

We  will  transmit  a  number  of  waveforms  from  a  speaker  and  measure  .he  received  waveforms  a.  a 
number  of  locations.  The  types  of  waveforms  that  we  will  use  include  the  followm-. 

1  Tonals,  i,e.,  waveforms  consisting  of  a  single  frequency.  We  expect. o  use  between  five 
and  ten  different  tonals  in  the  fiequency  range  from  20  Hz  to  1  KHz.  These  waveforms 
will  be  used  to  characterize  die  rime-varying  behavior  of  the  channel  both  in  terms  of 
amplitude  and  phase  fluctuations. 

a  Muld-tonals,  i.e.,  waveforms  consisting  of  multiple  frequencies.  We  will  combine  the 
single  tone  waveforms  to  generate  the  muld-tone  waveforms.  These  waveforms  wi 
be  used  to  characterize  the  time-varying  behavior  of  the  channel  and  how  fluctuations 
in  at  one  frequency  can  be  related  to  fluctuations  at  a  different  frequency. 

3  Pulse  trains.  These  waveforms  will  be  used  to  simulate  sources  of  potential  interest 
'  We  will  use  both  single  period  pulse  trains  and  multiple  period  pulse  trains.  The  funda¬ 
mental  frequency  of  the  pulse  trains  will  be  between  10  Hz  and  100  Hz. 

4.  Broadband  modulated  waveforms.  These  waveforms  will  allow  time-varying  mea 
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suremem  of  .he  channel  impulse  response.  We  will  use  modulated  signals  wid.  fre- 
quency  ranges  from  20  Hz  up  to  1000  Hz. 

Each  wavefonn  will  be  appreximately  a  minute  in  duradon  and  will  be  preceded  wid,  a  standard 

header  to  mark  it.  . 

Measuren^en.  wiU  be  made  a.  a  ""  ^^ievl  Sd"::  "will 
phone  arrays.  Most  measurements  ^  n,ade  morning,  noon. 

^d^etntoTe^ach^f  the  days  of  the  experiment.  Most 
a  selected  subset  of  the  sensors. 

r=“EHSHfaS=S5 

whh  one  minute,  five  minute,  and  twenty  four  hour  averages  dunng  the  expenment. 

3  Test  Plan 

Aet  r^hmi-nt  A  fASSCE  Test  Plan;  Acoustic  Self  Location)  is  the  formal  test  plan  for  this  data  col- 

vide  calibrated  waveforms,  rather  than  the  propane  cannons. 


Activity 

Start  Date 

End  Date 

Arrive 

Nov  30 

Nov  30 

Setup  Equipment 

Dec  1 

Dec  4 

Collect  Data 

Dec  7 

Dec  10 

Pack  Equipment 

Dec  11 

Dec  11 

IdDie  1.  Icai  piaii 

The  sneaker  for  the  wavefoims  will  be  placed  on  the  South  Tower  (page  4  of  Attachment  A)  tmd 

4  Planned  Analysis 

vses: 
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1.  Signal  to  noise  ratio  analysis  to  verify  data  quality. 

2.  Spectral  analysis  to  quantify  basic  signal  propagation  charactenstics. 

3.  Channel  impulse  response  extraction  to  get  a  quick  look  at  the  channel  charactenstics. 

Laboratory  analysis  will  focus  on  the  use  of  nonlinear  systems  modeling  techniques,  including: 

1.  Mutual  information  analysis  to  determine  the  proper  time  delay  for  subsequent  analy- 
sis. 

2  Time-delay  embedding  and  state  space  reconstmction.  These  techniques  are  used  to 
determine  if  the  data  can  be  modeled  as  a  low  dimensional  nonlinear  sy  stem.  Measure¬ 
ments  of  the  time-varying  amplitude  and  phase  response  of  the  channe  an  o  t  e 
time-varying  impulse  response  will  be  used. 

3.  Hidden  Markov  Models  as  a  method  for  signal  detection  in  fluctuating  channels.  We 
will  train  signal  detectors  based  on  observed  characteristics  and  test  these  detectors  to 
determine  their  ability  to  enhance  signal  detectability. 

Analyses  will  be  correlated  with  environmental  conditions  with  direction  of  arrival  (for  the  circu¬ 
lar  array)  and  with  distance  from  the  source. 
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Modificstions  to  ASSCE  Tost  Plsn 


1  S©nsor  Locations 

•  »t,  A  c  c  n?  tf»ct  nlan  call  for  two  circular  arrays  and  eight  ground  sensors. 
The  specifications  in  the  ASSCE  pl^ 

The  following  modifications  are  made  to  the  AbSCt  lesi  pian  rcg  5 
for  the  ARO  data  collection: 

.  The  circular  array  located  near  the  North  Tower  will  be  used. 

•  The  eight  ground  based  sensors  will  be  used. 

.  The  circular  array  located  near  the  South  Tower  wUl  be  used  but  is  only  of  second^ 
interest.  If  the  circular  array  at  the  Noith  Tower  fails  then  the  circular  array  from  t  e 
South  Tower  can  be  used  as  a  replacement.  . 

.  Thr-  micmphones  will  be  placed  on  the  North  Tower  at  heights  of  approximately  30, 
S,  and  90  teK.  These  microphones  will  be  mounted  inside  styroioam  to  reduce  wi 

noise. 

2  DAPS  Channel  Assignments 

During  the  recording  of  the  ARO  datasets  the  channel  assignments  on  the  Ncnh  Tower  DAPS 
will  be  changed  as  follows: 


Channel 

Group  1 
Sensor 

Group  2 
Sensor 

0 

8  foot  array  0 

remote  0  -  3  km 

1 

8  foot  array  1 

remote  1  -  2  km 

2 

8  foot  array  2 

remote  2  - 1 .5  km 

3 

8  foot  array  3 

remote  3  - 1  km 

4 

8  foot  array  4 

remote  4  -  500  ft 

5 

8  foot  array  5 

tower  0-30  feet 

6 

8  foot  array  6 

tower  1-60  feet 

7 

8  foot  array  7 

tower  2  -  90  feet 

laDie  I ;  . . 

This  chan.e  implies  that  channels  5, 6,  and  7  of  Group  2  will  be  changed  from  remote  5, 6,  and  7. 
Led  in  the  self-tocaliaation  tests  to  tower  microphones  used  in  the  propagation  tests. 


3  Test  Setup 

In  addition  to  the  Bst  setup  specified  in  the  ASSCE  plan  we  wiU  make  the  followuig  addrnons 
•  Day  2  -  test  ASL  speaker  connections  and  playback 
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.  Dav  3  (after  DAPS)  setup  -  test  propagation.  A  special  test  setup  tape  will  be 
She^rd^DipS  wavefoL  will  be  examined  in  the  field  and  FedEx  to  Loeltheed 

Sanders  for  analysis. 

4  Recording  Conditions 

The  default  recording  conditions  wUl  be  with  a  500  Hz  low  pass  filter  (eight  order  Butterw^) 

H  ramolina  rate  of  2083  Hz  (same  as  for  the  self-localization  expenments).  If  a  selectable  low 
Tass  r  is  Sav^Lwe  tale  for  the  test  then  alternative  selections  will  be  specified  accord- 

ing  to  the  capabilities  of  the  filters. 

5  List  of  Waveforms 

Fanh  waveform  to  be  transmined  will  consist  of  the  following  format: _ 


Description 

Center  Frequency  (Hz) 

Duration  (secc 

13  chip  Barker  code 

100 

0.1 

13  chip  Barker  code 

200 

0.1 

13  chip  Barker  code 

300 

0.1 

1 3  chip  Barker  code 

400 

0.1 

Silence 

0 

0.5 

13  chip  Barker  code 

250 

0.1 

Five  bit  waveform  id 

250 

0.5 

13  chip  Barker  code 

250 

0.1 

Silence 

0 

1.0 

Data 

variable 

variable 

Silence 

0 

10.0 

Table  2:  Basic  waveform  formating. 

The  first  four  Barker  codes  are  used  to  generate  a  timing  mark.  The  total  bandwidth  of  the  signal 
used  for  the  timin-  mark  is  approximately  500  Hz  so  the  timing  resolution  will  be  approximately 
S  0%  Londs  Th"e  waveform  id  sequence  will  consist  of  five,  0.1  second  duration  segments  each 
of“  m  eilr  contain  a  13  chil,  Barker  code  or  silence.  This  id  capabiiity  wrll  provrie  or  up 
to  32  distinct  waveforms.  The  total  duration  of  the  header  and  the  silence  segments  is  -. 

ends.  The  data  sequences  are  as  follows:  _ 

r;^ - ~  ^  Description  I  (seconds) 

Tl  Sii^e  Silence  I  ^0 - - 

~T~  Tone _ _ _ _ ^ - - 

Table  3:  Waveform  descriptions. 


Duration  (seconds) 
60 
60 
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Tone 

Tone 

Tone 

Tone 

Tone 

Tone 

Tone 


Description 


40  Hz 
60  Hz 
80  Hz 
100  Hz 
150  Hz 
200  Hz 
300  Hz 


Lockheed  Sanders,  Inc. 


Duration  (seconds) 


9 

10 

Tone 

500  Hz 

11 

Multi-tone 

20  Hz  and  40  Hz 

12 

Multi-tone 

20  Hz  and  80  Hz 

13 

Multi-tone 

40  Hz  and  80  Hz 

14 

Multi-tone 

100  Hz  and  150  Hz 

15 

Multi-tone 

1 00  Hz  and  200  Hz 

16 

Multi-tone 

100  Hz  and  300  Hz 

17 

Multi-tone 

100  Hz  and  400  Hz 

18 

Muiti-tone 

1 00  Hz  and  500  Hz 

19 

Pulse  Train 

Fundamental  =  15  Hz 

20 

Pulse  Train 

Fundamental  =  30  Hz 

21 

Pulse  Train 

Fundamental  =  65  Hz 

22 

Multi-Puise  Train 

Fundamentals  =  15  Hz  and  65  Hz 

.  •  .  1  I-,  CC  Ul'T  1 

23  Mulll-Pulse  Train  Fundamentals  -  au  m  a.. u  _ _ _ 

17~0PSK  Modulation  Baud  Rale  -  50  Hz.  Center  Freque^  ■  250  _ 60_ - - 

IT  QPSK  Modulation  Baud  Rale  -  100  Hz,  Center  Frequency  ■  250  Hz _ 60 - 

26  QPSK  Modulation  Baud  Rate  -  250  Hz,  Center  Frequency  -  250  Hz - - - 

27  QPSK  Modulation  Baud  Rate  =  50  Hz,  Center  Frequency  -  1 00  Hz - 

28  I  QPSK  Modulation  Baud  Rale  .  100  Hz,  Center  Frequency  -  100  Hz_ - 60 - 

iriilPSK  Modulation  |  Baud  Rale  ■  50  Hz,  Center  Frequenq^OO  Hz _ 60 - 

an  I  mpdK  Mod, nation  I  Baud  Rate,  too  Hz,  Center  Fiequenc^40^ - 60 - 1 

- - -  Table  3:  Waveform  descriptions. 

AU  signals  will  be  generated  at  a  2083,33  Hz  sampling  rate  and  low  pass  filtemd  with  a  800  Hz 


tesl_plan.fm 


6 


November  30, 1S92  6:44  pm 


ARO  Test  Plan 


Lockheed  Sanders,  Inc. 


low  pass  filter  prior  to  digitization. 


Table  4:  Data  collection  play  list. 


The  play  list  for  the  special  test  tape  will  be  as  follows: 


Segment  Number 

Waveform  ID 

1 

0 

Table  5:  Test  tape  play  list. 
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Abstract 


This  plan  describes  the  objectives,  procedures,  equipment,  schedule  aud  level 
of  effort  needed  to  gather  experimental  data  for  investigating  thafwm 

and  self-orientation  using  propane  cannon  sources.  ^ 

be  performed  concurrently,  is  the  ARO  experiment  which  designed  to  collect 
propagation  data  using  speaker  sources.  The  purpose  of  the  A  JL  i 

genLate  practical  models  for  the  spatially  and  temporally  varying  acoustic  ch^ 
fn  the  deLt  environment.  Raw  data  wUl  be  recorded  in  a  field  test  environment 
near  the  DIRT  Site  at  White  Sands  Missile  Range.  This  data  will  be  used  o 
demonstrate  self-location  algorithms  in  the  Lockheed-Sanders  SPCOT  laboratory 

at  a  later  date. 


1 


1  Objective 

The  fundamental  objective  of  this  exercise,  comprised  of  a  field  test,  is  to  gather  data 
to  be  used  to  demonstrate  self-location  algorithms  which  demonstrate  the  ability  to 
eliminate  the  requirement  for  beamformer  calibration,  self-righting,  north-seeking,  ^d 
position  location  capabilities  in  each  node  of  a  distributed  acoustic  sensor  system.  The 
development  of  such  a  technique  could  impact  systems  such  as  WAM,  netted  ^nes,  and 
distributed  acoustic  sensors,  by  reducing  the  size,  weight,  and  cost  of  each  node  through 
the  elimination  of  self-righting  capability,  compasses,  and  position  electronics. 

This  self-location  objective  will  be  met  by  the  following:  Collection  of  data  during  a 
field  experiment  from  both  circular  sensor  arrays  and  remotely  located  sensors,  inclu  ing 
ambient  noise  levels  and  impulsive  blasts  with  at  least  four  impulsive  cannons  placed  at 

different  locations  in  the  test  area.  •  r  j  * 

The  objective  of  the  ARC  experiment  will  be  met  by  the  following:  Collection  of  data 
at  various  times  of  the  day  that  will  consist  of  synthetic  signals  including  harmonically 
related  tones  and  broadband  signals.  These  signals  will  be  broadcast  from  the  Servo- 
drive  speaker/amplifier  system  that  is  mounted  in  the  South  Tower.  The  signals  will  be 
broadcast  for  no  longer  than  15  minute  durations  at  intervals  in  the  morning,  midmornmg, 

noon,  and  midafternoon.  ^i,  j  + 

Whenever  possible,  a  Sun  workstation  will  be  used  to  perform  analysis  on  the  data 

and  also  to  exercise  the  demo  algorithms. 


2  Approach 

The  data  collection  exercises  will  be  performed  outside  in  a  field  test  situation.  The 
test  will  be  performed  at  the  DIRT  Site  at  the  White  Sands  Missile  Range,  near  the 
Orogrande  Range  Camp. 


2.1  Field  Test 

The  test  will  consist  of  recording  a  series  of  signature  measurements  made  up  of  ambient 
noise  and  impulsive  noise  from  cannons.  Up  to  five  single  microphones  (remote  sensors) 
will  be  located  along  a  straight  line  north  and  south  of  the  DAPS  located  at  the  North 
tower.  One  of  these  single  sensors  will  be  placed  near  the  circular  array  at  the  North 
tower  site.  Additionally,  three  remote  sensors  will  also  be  located  along  a  parallel  line 
2km  to  the  west.  All  the  sensors  will  be  connected  to  the  DAPS  via  field  wire.  This 
field  wire  will  be  supplied  by  ARDEC.  Additionally,  ARDEC  wiU  provide-  a  field  wire 
laying  apparatus  that  is  man  portable  and/ or  a  field  wire  laying  kit  that  attaches  to  the 
HMMWV.  ARDEC  will  also  provide  a  kit  that  enables  the  field  wire  to  be  rolled  up  at 
the  conclusion  of  the  test,  if  required.  ASL  wil  provide  a  HMMWV  and  will  investigate 
any  permissions  required  to  drive  the  vehicle  both  on  and  off  road  in  the  general  area. 


2 


A  speaker  »ill  be  placed  in  the  South  tower  for  the  ARO  experiment,  A  second  DAPS 
with  a  circular  array  will  also  be  located  at  the  South  tower 

With  respect  to  the  ARO  experiment,  an  acoustic  source  (lou  p  )  a- 

on  the  South  MET  tower  approximately  100  feet  above  ground.  A  senes  acoustic 
sensors  will  be  placed  on  the  ground  at  var^ng  distances 

speaker:  a  single  sensor  at^O  ^^-^"uTthe  !eto!"s  fetW^part  of  tht 

expeJtoent  'vriU  be  connected  io  the  DAPS  at  the  North  tower.  Synthetic  signals  mcludmg 

h^monically  related  tones  and  broadband  signab  will  be  Wdcast  from  the 

r^iei  at  each  of  the  sen^r  locations.  These  data  wffl  then  be  verified  on  a  portable 

workstation  in  the  field. 

2.1.1  Array  Data  Collection 

The  equipment  for  this  part  of  the  test  will  be  arranged  as  in  Figure  1.  Each  DAPS  will 
support  Ip  to  two  groups,  each  consisting  of  eight  channels  “P  “ 

betLn  primary  and  remote  microphone  sensors.  The  D^  DAPS  ^ 
the  North  tower  and  will  support  a  circular  array  and  eight  remo  e  array  . 

DAPS  will  be  located  at  the  South  tower  and  wiU  support  the  circular  array  In  orde 
io  make  future  processing  of  data  from  a  given  DAPS  easier,  only  sensors  of  one  type 
(either  high  or  low  gain)  will  be  used  in  the  test. 

2.1.2  Test  Conduct 

The  test  will  be  conducted  at  the  White  Sands  Missile  Range  NM  during  the  time 
period  timber  30,  1992  to  December  12,  1992.  The  test  is  eduled  to 
Lys  The  first  day  will  be  for  equipment  transportation  to  the  field  test 
freight  receiving  area.  The  next  four  days  will  be  for  equipment  setup  and  checkout.  Th 
MowlnTSur  days  will  be  for  data  collection,  and  the  last  ‘ 

tear-down  and  transportation  back  to  the  shipping  oca  ion.  11+  ^  data  will  be 

be  conducted  at  various  times  throughout  the  day  as  needed.  The  collected  data  will 

checked  for  data  integrity  during  the  test 

data  will  be  sorted  and  annotated  to  maximise  future  usabihty.  The  following  sections 
give  detailed  information  on  each  of  the  stages  of  the  field  test. 

Test  Equipment  Transportation/Setup  This  portion  of  the  field  test  will  ac- 

complished  in  five  days.  Day  one  will  consist  of  ^  ^ 

test  site  A  U-Haul  truck  or  another  appropriate  vehicle  will  be  rrated  for  a  day  to 
mLport  the  material  from  the  air  freight  terminal  in  El  Paso  to  the  t«t  range  The 
day  will  be  spent  unpacking  the  equipment  from  the  shipping  containers,  loaAng  it  on  o 
thi  vehicle  driving  out  to  White  Sands,  NM,  and  then  unloading  the  equipment  into 
ifuipmeWpefsonnel  trailer.  A  fork  lift  to  unload  the  boxes  from  the  truck  will 
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•  8  foot  array 
■  Tower 

#  Remote  Sensor 


Figure  1:  General  Site  Setup 


be  required  during  setup.  One  wiU  also  be  required  during  teardown  and  shipment  of 

material  back  to  the  freight  termiiial.  ,  ,  i  .  . i  • 

The  Mowing  four  days  will  be  used  to  setup  and  checkout  the  equipment  at  the  field 

test  site.  The  workstation  wiU  be  set  up  in  the  equipment  trails.  Also,  the  circular 
arrays  which  wiU  be  located  at  the  North  and  South  Towss  will  be  cons  lurted^  A 
signfficant  task  in  terms  of  time  will  be  that  of  laymg  out  the  remote  sensor  field  wir  _ 
Trends  of  the  field  wire  will  be  finished  with  XLR  three  conductor  type  “minors  (^if 
the  connectors  have  not  already  been  attached  -  plug  on  one  end  and  a  so<±et  on  th 
other  end).  The  sensors  that  will  be  placed  in  the  path  of  the  speaker  P«>P^8^  ‘°“;Lnce 
the  North/South  Tows  line)  will  be  laid  out  neict  to  the  <hrt  road  with  the  ^tsi  ance 
of  a  vehide.  The  remote  sensor  cables  on  the  line  located  2km  away  and  paraUel  to 
the  Nortb/South  Tows  Hue  will  also  be  laid  out  with  the  assistance  of  a  cross-country 
vehicle,  if  one  is  available.  Trenche(s)  wUl  be  dug  in  the  dirt  road  where  necessary  to 

prevent  damage  to  the  field  wire  due  to  road  traffic. 

Also,  during  the  setup  stage  both  the  remote  and  circular  ^ay  sensor  locations  will 

be  surveyed  and  recorded.  This  will  be  performed  with  the  aid  of  a 
finder  and  also  a  differential  GPS.  The  speaker  and  amplifier  system  will  be  mounted  up 
in  the  tower  bv  ASL,  along  with  the  remote  sensor  that  is  to  be  positioned  m  the  tower. 
Once  setup  is  Completed,  a  system  check-out  wiU  be  done,  including  verification  of  remote 
sensors,  local  sensors,  D.APS,  Sun  workstation,  speaker  system,  and  propane  cannons. 
See  the  Appendix  for  a  more  detailed  schedule  of  the  test  setup  and  data  collection 
The  DAPS  units  will  be  stored  in  the  test  trailer  at  night  and  the  arrays  wi  e 

wrapped  in  weatherproof  tarps  and  stored  underneath  the  trailer.  The  remote  sensors 
will  be  also  be  housed  in  the  trailer  at  night  as  well,  however,  the  field  wire  cabhng  wi 
be  left  out  for  the  duration  of  the  experiment. 


Recording  Assignments  Tables  1  and  2  give  sensor  position/DAPS  recording  as- 
signments.  These  recording  assignments  may  change  due  to  circumstances  of  the  test  as 

required. 

Ambient  Noise  Data  Collection  This  portion  of  the  test  is  designed  to  record  a  series 
of  noise  level  measurements  of  the  environment.  The  DAPS  will  be  used  as  the  recor  ng 
device.  The  data  will  be  collected  to  provide  a  baseline  to  verify  proper  operation  of  la  er 
portions  of  the  experiment.  Ambient  noise  measurements  will  be  recorded  throughout 
the  test  to  provide  data  under  various  propagation  conditions. 

Impulsive  Event  Data  Collection  This  portion  of  the  test  is  designed  to  record  a 
series  of  impulsive  events  from  the  orchard  cannons.  These  cannons  wiU  be  placed  at 
locations  within  the  test  site  and  their  locations  wiU  be  determined  via  the  l^er  range¬ 
finder  and  differential  GPS.  Acoustic  recordings  will  be  made  with  all  of  the  cannons 
firing  at  a  known  average  rate,  but  not  synchronized  in  any  way.  This  portion  of  the 
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Table  3:  Test /Schedule  Cross-Reference 


Test  1 

Basehne  AM 

Test  2 

Basehne  PM 

Test  3 

Angular  Resolution  AM 

Test  4 

Angular  Resolution  FM 

Test  5 

Ranging  Resolution  AM 

Test  6 

Ranging  Resolution  PM 

Test  7 

Small  Constellation  AM 

Tests 

Small  Constellation  PM 

Test  A 

TBD 

Test  B 

TBD 

experiment  wiU  take  approximately  four  days.  Daily  maintenance  wiU  be  performed  on 

the  remote  power  supplies  as  is  needed.  . 

A  number  of  different  configurations  will  be  used  as  data  c(^ection  scenarios. 
wiU  be  coUected  for  as  many  of  these  scenarios  as  is  possible.  The  initial  configura  ion 
will  serve  as  a  baseline  configuration  and  wiU  utilize  as  many  of  the  sensors  as  possmle. 
Data  will  be  coUected  in  the  baseline  configuration  over  a  several  hour  period  including 
morning  and  evening  if  possible,  to  provide  comparisions  on  performance  under  various 
climactic  conditions.  Another  configuration  wiU  be  set  up  to  test  the  resolving  capabi  ities 
of  the  axrav  for  varied  spacing  between  the  sources.  The  sensor  arrangement  will  be 
similar  to  the  baseline  test,  but  the  sources  wiU  be  placed  at  increasing  distances.  A 
similar  test  wiU  be  run  to  investigate  ranging  resolution  of  the  system.  The  sources  wi 
be  placed  as  increasing  distances  in  a  linear  fashion  outside  of  the  sensor  area  and  data 
wiU  be  gathered.  K  time  permits,  the  non-colUnear  sensors  wiU  be  moved  closer  topther 
to  investigate  the  effects  of  a  closer  constellation  of  sensor  elements.  For  this  expenment, 
the  source  arrangement  from  the  angular  resolution  test  will  be  utiUpd. 

The  baseline  sensor /cannon  arrangement  is  shown  in  Figure  2.  The  sensor/ cannon 
angular  resolution  test  arrangement  is  shown  in  Figure  3.  The  sensor/cannon  setup  for 
the  ranging  configuration  is  shown  in  Figure  4,  and  the  setup  for  the  small  constellation 
is  shown  in  Figure  5.  Table  3  shows  a  cross-reference  between  the  test  configuration  and 
the  test  scedule.  Note  that  optional  tests  A  and  B  will  be  determined  dunng  the  test  and 
wiU  either  be  used  to  collect  data  in  one  of  the  previously  mentioned  configurations  or 
wiU  be  used  to  collect  data  in  a  configuration  which  may  be  desirable  due  to  particular 
topographical  irregularities  at  the  test  site. 


ARO  Data  Collection  The  ARO  experiment,  which  will  be  performed  concurrently, 
is  designed  to  collect  propagation  data  using  a  speaker  source.  The  purpose  of  this  exper¬ 
iment  is  to  generate  practical  models  for  the  spatially  and  temporaUy  varying  acoustic 
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Figure  2:  Baseline  Configuration 
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Figure  3:  Angular  Resolution  Configuration 
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Figure  4:  Ranging  Resolution  Configuration 
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Figure  5:  Small  Constellation  Configuration 


channel  in  the  desert  environment. 

Data  will  be  collected  at  various  times  of  the  day  and  will  consist  of  synthetic  sig¬ 
nals  including  harmonically  related  tones  and  broadband  signals.  These  sipials  wil  e 
broadcast  from  the  Servo-drive  speaker/amplifier  system  that  is  mounted  in  the  South 
Tower.  The  signals  will  be  broadcast  for  no  longer  than  15  minute  durations  at  intervals 
in  the  morning,  midmorning,  noon,  and  midafternoon.  Time  permitting,  these  scenarios 
will  be  recorded  over  the  course  of  more  than  one  day  of  the  test.  The  propane  cannons 

will  need  to  be  turned  off/on  during  these  propagation  tests. 

This  experiment  will  run  for  the  duration  of  the  self-location  test,  and  the  time  of  day 
when  the  synthetic  signals  will  be  broadcast  may  be  modified  as  appropriate  to  coincide 
with  down-time  of  the  other  experiment.  The  sensor  configuration  to  be  used  will  be  the 
same  as  the  baseline  sensor  configuration  used  in  the  impulsive  data  collection  exercises. 

Test  Cleanup  This  portion  of  the  test  will  take  approximately  one  day.  The  DAPS 
equipment,  circular  arrays,  and  remote  sensors  will  be  gathered  up  and  returned  to  the 
site  trailer.  The  equipment  will  then  be  boxed  up  and  returned  to  the  shipping  site  in 
El  Paso.  Again,  a  U-Haul  type  vehicle  will  be  rented  to  transport  the  material  to  the  air 

freight  terminal. 

2.1.3  Test  Personnel 

The  following  personnel  are  required  to  run  the  test.  A  brief  description  of  their  duties 
is  listed  below. 

Test  Coordinator  Responsible  for  the  overall  implementation  of  this  test  plan  and 
coordination  of  the  various  aspects  of  this  test.  Makes  decisions  as  required  for  running 
the  tests.  Will  also  assist  in  data  collection  exercises  as  well  as  set-up  and  tear  down. 

DAPS  Equipment  Operator/Scribe  -  2  Responsible  for  the  operation  of  the  DAPS 
unit.  Ensures  that  the  unit  is  operating  properly  prior  to  start  of  the  test,  and  maintains  it 
during  the  test.  Records  the  configuration  of  the  system  at  test  time  (channels,  sampling 
rates,  etc).  Makes  sure  that  recording  tapes  are  properly  loaded  and  labeled  Mak^ 
sure  system  clock  is  properly  synchronized  and  set.  Records  information  in  the  log  book 

including  the  following: 

(a)  Impulse  Times  -  Records  the  time  that  the  impulsive  events  were  heard  at  the 
DAPS. 

(b)  Interference  -  Records  when  any  interference  (i.e.  vehicles,  planes,  helos,  etc.)  are 
present  in  the  test  area. 

(c)  Start /Stop  Times  -  Record  the  start  and  stop  times  of  each  portion  of  the  experi¬ 
ment  and  any  pertinent  details  associated  with  each  data  collection  phase. 


(d)  Location  -  Record  the  location  of  all  remote  and  local  sensors. 


(e)  Other  -  Record  any  other  pertinent  information. 


wy  - - 

All  personnel  at  the  site  will  assist  in  setup  and  maintenance  of  the  data  collection  and 
signal  source  equipment.  This  includes  daily  setup  of  the  propane  cannons,  microp  on 
arrays  recording  units,  remote  microphone  power  supplies,  etc.  As  time  permits,  the  tes 
personnel  will  perform  analysis  on  the  data  that  has  been  collected.  This 
will  be  performed  on  the  Sun  workstation  and  will  be  used  to  provide  feedback  as 

proper  system  operation  and  configuration. 


3  Description  of  Equipment 

3.1  DAPS 

The  main  instrumentation  recorder  to  be  used  in  this  test  is  Sanders’ Digital  Acquisition 
and  Processing  System  (DAPS)  developed  by  the  NCTR  lab.  Two  DAPS  units  will  b 
used  a  24Vdc  “desert  DAPS”  configuration  and  a  DAPS  m  the  TDOA  con  gura  ion. 
In  the  field  test,  the  Desert  Daps  will  be  connected  to  the 

element  circular  array  located  at  the  North  tower.  The  TDOA  DAPS  will  be  conne 
to  an  eight  element  circular  array  located  at  the  South  tower.  r  _ 

Two  analog  modules  (housed  in  the  VME  chassis)  provide  the  anti-ah^ing  filters 
each  channel  prior  to  digitization  by  the  DAPS.  The  DAPS  is  controlled  by  a  Motorola 
68020  CPU  and  is  able  to  store  up  to  2  Gbytes  of  signal  information  at  a  rate  of  20 
Kbytes/second.  The  storage  medium  is  8  mm  video  cassette  tape. 

The  DAPS  records  a  time  stamp  on  the  data  from  its  own  time  o  _ 

the  microphones  are  on  different  groups  but  are  recorded  with  the  same  DAPS  a  0^5 
second  ambiguity  in  synchronization  may  exist.  Written  data  will  provide 
the  tests.  An  oscilloscope  may  be  used  to  check  the  integrity  of  the  sensor  data  in  the 

field. 


3.2  Signal  Sources 

For  generation  of  the  impulsive  noise  signals,  up  to  five  a^rtillery  noise-makers  wUl  be 
used.  Four  of  these  systems  will  be  provided  by  the  Signal  Processing  en  ° 
nology,  Lockheed-Sanders.  ASL  will  provide  units  as  necessary  to  function  as  additional 
impulsive  noise  sources  and  also  to  serve  as  replacement  units  m  the  event  that  spares 

For  generation  of  the  synthetic  harmonically  related  tones  and  broadband  signa  s, 
a  high  power  speaker  system  will  be  used  in  conjunction  with  a  tape  deck.  ASL  wiU 
supply  the  Servo-drive  speaker  and  the  high  power  amplifier  while  Lockheed-S^d«s 
wifi  provide  the  tape  deck  and  pre-recorded  synthetic  signals.  ASL  mil  provide  the 
harness  for  mounting  the  speaker  in  the  tower  as  well  as  assistance  mounting  the  speaker 


11 


system  in  the  South  Tower.  Lockheed-Sanders  will  need  access  to  the  speaker  system  at 
various  times  throughout  the  test. 

3.3  Radio  Communications 

Lockheed-Sanders  will  bring  two  hand-held  radios  to  f 

cations  both  during  site  setup  and  during  test  exercises  to  12/11/92)  These 

radios  operate  at  158.280MHz  and  153.395MHz  and  are  Licensed  by  the  FCC  opera¬ 
tion  in  the  US.  The  emission  designator  for  these  radios  is  20K0F3E.  The  call  le  ters 
these  radios  are  WNMU754  and  WRP698,  respectively. 

3.4  Surveying  Tools 

ASL  will  provide  a  Laser  Transit  for  surveying  the  sensor  locations  and  the  cannon 
locations  ASL  will  also  provide  assistance  in  surveying  these  points  during  setup  as  wel 
Z  intermadlL  survey  points  required,  if  the  Laser  Trartsifs  effect.ve  survey 
range  is  shorter  than  the  distance  from  the  reference  point  to  the  sensor  and/or  cannon 
locations.  If  additional  points  need  to  be  surveyed  during  the  data  collection  portion  o 
the  experiment,  either  ASL  personnel  or  Lockheed-Sanders  personnel  will  perform  the 
site  survey  with  the  Laser  Transit.  Up  to  twelve  cannon  locations  and  twelve  sensor 

location  will  need  to  be  surveyed.  .  ,  •  j 

Lockheed-Sanders  will  provide  2  portable  GPS  systems  to  survey  points  in  a  pseudo¬ 
differential  GPS  mode.  This  surveying  will  be  done  independently  from  the  surveying 
that  is  performed  with  the  Laser  Transit  and  is  intended  to  verify  the  survey  points  of 

the  Laser  Transit. 
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Appendix  A:  Equipment  list 

(1)  24Vdc  Digital  Acquisition  and  Processing  System 

(1)  24Vdc  Digital  Acquisition  and  Processing  System  (TDOA) 

(1)  Power  Supply  for  TDOA 

(4)  12  Volt  Deep  Cycle  Marine  Batteries 

(2)  Portable  GPS 

(4)  M-4  Single  Detonation  Cannon 
(1)  RF  Modem 

(1)  Set  Battery  Connectors 

(2)  Exabyte  Tape  Drive  in  External  Shoebox 

(1)  Exabyte  Tape  Drive  (spare) 

(2)  Random  Console  terminal  with  charging  cable 

(1)  Compass 

(1)  Sensor  Calibration  Kit 

(1)  Tape  Deck  to  play  synthetic  signals  ^  c  r 

(2)  Servodrive  speaker  and  associated  amplifier  (Gov’t.  Supplied) 

(1)  Speaker  Harness  for  tower  mounting 

(1)  South  Tower 

(30)  Countryman  Microphones 
(25)  Single  Sensor  Windscreens 

(2)  l’x.5’x.25’  styrafoam  squares  (wind  screen  for  tower  mounted  microphone) 
(2)  2’  Velcro  mounting  straps  (for  tower  mounted  microphone) 

(2)  3”  (approx)  pipe  clamps  (for  tower  mounted  microphone) 

(1)  120’  DAPS-microphone  cable  (for  tower  mounted  microphone) 

(2)  8’  Circular  Arrays  and  construction  hardware 
(4)  Junction  Boxes 

(8)  Remote  Sensor  Phantom-Power  Supplies 

(150)  9  Volt  Batteries  for  Remote  Sensor  Power  Supplies 

(3)  2km  Rolls  of  Field  Wire 
(2)  1km  Rolls  of  Field  Wire 

(1)  1.5km  Roll  of  Field  Wire  (ARO) 

(1)  2km  Roll  of  Field  Wire  (ARO) 

(1)  3km  Roll  of  Field  Wire  (ARO) 

(2)  0.5km  Rolls  of  Field  Wire 
(1)  150’  Roll  of  Field  Wire 
(1)  Field  Wire  Repair  Kit 
(1)  Cable  Laying  Apparatus 

(3)  Boxes  8mm  Data  Tapes 

(1)  Battery  Recharger 
(12)  Propane  Tanks 

(2)  Pairs  earmuffs 

(1)  Adjustable  Ladder 
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(5)  Weatherproof  Tarpaulins 

(1)  SparcStation  SLC  with  8inm  tape  drive  and  external  hard  disk 

(2)  P-100  Hand-Held  Radios  and  rechargers 
(1)  Oscilloscope 

Duct  tape,  elastic  cords,  cabling,  spikes,  mallet 
Dust  Shields  for  DAPS  and  Exabyte  Tape  Drive 
(1)  Electroline  Duct  Seal 

Assorted  small  tools  (screwdrivers,  etc),  power  cords,  etc 
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Appendix  B:  Test  Schedule 
Week  1 


Day  1 
Day  2 


Day  3 


Day  4 


1/2  day  -  transport  material  out  to  test  site,  begin  unpacking 
1  /4  day  -  quick  site  survey 

3/4  day  *  begin  laying  field  wire,  attaching  connectors 

1  /2  day  -  laying  remaining  field  wire 
1  /2  day  -  begin  surveying  sensor  locations 

entire  day  to  survey  remaining  sensor  locations 
and  start  surveying  cannon  locations 


Day  5 


1/2  day  -  survey  remaining  cannon  locations 
1/2  day  -  debug/repair  as  necessary 


Week  2 


Day  1 
Day  2 
Day  3 
Day  4 
Day  5 

Day  6 


1  /2  day  -  Test  1 
1  /2  day  -  Test  2 

1  /2  day  -  Test  3 
1  /2  day  -  Test  4 

1  /2  day  -  Test  5 
1  /2  day  -  Test  6 

1  /2  day  -  Test  7 
1/2  day  -  Test  8 

1  /2  day  -  Optional  Test  A 
1/2  day  -  Optional  Test  B 
OR  .  .  , 

entire  day  to  pack  and  ship  material  back  to 
Air  Freight  terminal 

entire  day  to  pack  and  ship  material  back  to 
Air  Freight  terminal  if  not  done  prior  day 
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Final  Report 


SPCOT 


MEMO:  August  25, 1993 


Edward  C.  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.Real 

885-5147 

RE: 

Elimination  of  interfering  signals  from  acoustic  data. 

FAX  885-0631 

CC: 

I 

I 

I 

I 


TMs  is  a  summary  report  on  my  efforts  to  date  to  remove  the  interfering  signals  from  the 
ARO  data  set  The  basic  operarional  parameters  for  this  work  are  outlined  m  the  followmg 


Parameter 

Value 

Comments 

Sample  Rate 

2083.33  Hz 

As  given  in  the  ARO  test  plan. 

Data  Acquisi¬ 
tion  Date 

Dec.  7. 1992 

As  given  on  data  tape  label. 

lime 

1  to  1:10  pm 

ibid. 

Location 

North  Tower, 

WSMR 

ibid. 

Tape  #  That 
Data  Is  On 

2 

ibid. 

Channel  Used 

3 

Relatively  clean  data. 

Group  #  Used 

1 

There  are  2  groups  on  this  tape. 

.•  _■  ¥> _ 

Table  1:  Basic  Operational  Parameters 


For  this  study  the  foUowing  signal  types  were  extracted  by  hand  from  the  data  file 


identi 


Signal 

ID 

Class 

Description 

0 

Silence 

No  intentional  signals  emitted. 

18 

Multi-tone 

100  and  500  Hz. 

19 

Impulse  Train 

Impulses  emitted  at  a  15  Hz  rate. 

20 

Impulse  Train 

Impulses  emitted  at  a  30  Hz  rate. 

21 

Impulse  Train 

Impulses  emitted  at  a  65  Hz  rate. 

22 

Multi-Impulse 

Train 

Impulses  emitted  at  15  and  65  Hz  rates. 

23 

Multi-Impulse 

Train 

Impulses  emitted  at  30  and  65  Hz  rates. 

24 

QPSK 

Center  frequency  at  250  Hz. 

Baud  Rate  is  50  Hz. 

25 

QPSK 

Center  frequency  at  250  Hz. 

Baud  Rate  is  100  Hz. 

.  «T  1  X  .  Ol. _ 

Table  2:  Signal  Types  Used  In  This  Study 
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Preliminaries  -  The  Sounds  Of  Silence. 

The  goal  of  this  study  is  to  remove  unwanted  signals  from  the  enviroi^nt  without  1)  dis- 
tobig  the  signals  of  interest  and/or  2)  changing  the  chmacter  of  ^dato  m  s«oh  a  was^ 
to  distort  any  subsequent  measurements  of  the  acoustic  channel.  For  ^he  mos  p  , 

are  tones  generated  by  nearby  machinery;  althou^  this  is  not  necessa^- 
Uy  the  (only)  form  of  the  interference.  An  example  of  the  tonal  mterfermg  signals  is 
shown  in  figure  1. 

Fioure  1  is  a  time/ftequency  (spectrogram)  plot  of  a  bracketed  segment  of  silence  taten 
from  the  beginning  of  the  data  file  outlined  m  table  1.  The  tune  ongm  is  relahve  to  ^e 
beginning  offeplot,  and  intensity  (magnitude  in  dB)  is  identified  by  grey  scale.  Darter 
rejons  indicate  greater  arnpUmde.  Several  features  of  this  sUence  segment  are  readily  vis¬ 
ible.  and  warrant  further  discussion. 

The  first  set  of  features  we  consider  are  the  broad  band  signals  with  a  bandwidth  of  about 
4(X)  Hz  ranging  between  approximately  200  to  600  Hz.  These  are  evidenm  the  time  seg¬ 
ments  between  7  and  10  seconds,  and  again  between  80  and  83  seconds.  These  sisals  are 
riming  marks  which  were  emitted  to  bracket  the  silence  segment  which  lies  m  between 

(-10  to  80  seconds). 

Another  prominent  feature  is  the  300  Hz  tone  evident  from  about  83  to  90  seconds.  This  is 
the  tvmnnlng  of  sigual  8  which  immediately  foUows  signal  0  on  the  play  list. 

These  features  aside,  the  rest  of  figure  one  should  represent  sUence.  Qearly.  to  is  not  the 
case.  There  remain  several  distinct  tonal  artifacts  from  nearby  machinery.  There  tonah 
appear  to  be  harmonically  related,  with  a  frequency  spacmg  of  a  httle  over  100  Hz  A 
properly  centered  comb  filter  with  to  spacing  and  an  appropnate  set  of  gam  values  for 
the^th  might  suffice  to  eliminate  these  tones,  but  other  artifacts  m  die  noise  wodd  then 
require  separate  processing;  assuming  they  could  be  identified.  There  is  also  to  addiuonal 
problem  of  new  interferers  whose  spectral  characteristics  vary  over  time.  For  these  rea- 

Signal  0  :  Unprocessed  Silence 


30  40  50  60 

Time  (Seconds) 


Figure  1:  Spectrogram  of  unprocessed  silence. 
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sons,  it  is  mote  convenient  to  have  a  technique  which  ad^ts  to  the  data  at  hand  than  a  set 
of  rigid  filters. 

A  Projected  Solution 

The  three  major  requirements  on  the  interference  removal  techmque  are  then  1)  minimum 
impact  on  the  signals  of  interest,  2)  minimum  distortion  of  any  acoustic  channel  measure¬ 
ments  and  3)  data  adaptability.  These  three  requirements  by  no  means  uniquely  identify  a 
particular  tedmique  to  use.  There  are  many  teduuques  which  could  rnake  a  claiin  at 
achieving  these  goals.  However,  one  technique  which  recommends  itself  is  the  techmque 
of  null  space  projection. 

We  begin  by  gcgitning  that  the  interferers  are  stationary  over  a  long  enough  period  of  time 
that  a  cAgmpnt  of  the  data  containing  the  interferers  alone  may  be  characterized  via  a  pro¬ 
jection  operator  constructed  out  of  the  data  itself.  That  is,  let  be  a  segment  of  the  dam 
contrining  only  tb^  interfering  signal(s).  Obtain  a  set  of  n  such  data  segments  { (possi¬ 
bly  overlapping)  and  use  these  segments  as  a  basis  set  to  characterize  the  interferers.  The 
^cciimptifin  is  that  these  data  segments  {jc  }  are  linearly  independent,  so  it  is  important 
that  overlapping  data  segments  preserve  this  assumption.  This  requirenient  is  not  hard  to 
fulfill  in  practice.  Any  other  data  segment  (vector)  Xj  may  be  projected  into  this  space  via 
the  projection  operator 


P  A  [A"a]~^A^ 


(1) 


where  H  indicates  complex  conjugate  transpose  (Hermitian).  and  the  matrix  A  is  con¬ 
structed  from  then  element  data  set  {.r,}  as. 


A 


I  I  I 

.V,  x,  ...  x„ 

I  I  I 


(2) 


where  the  columns  of  A  are  made  up  of  elements  of  the  data  set  {.r,}  .  The  projected  data 
vector  X).  is  given  by 


A  refinement  to  this  technique  is  to  apply  a  reduced  rank  approximation  to  (2)  in  order  to 
select  only  those  components  of  A  which  contribute  the  most  to  the  projection.  The  ratio¬ 
nale  behind  this  approach  is  that  1)  the  strongest  contributors  to  the  projection  operator  are 
the  interferers  and  2)  we  want  to  distort  the  data  relating  to  the  acoustic  channel  as  little  as 
possible,  so  we  only  identify  that  which  we  don’t  want.  With  this  goal  we  reformulate  (1) 
in  terms  of  its  singular  vectors. 

The  Moore-Penrose  inverse  of  A  is  given  by 

/I*  -  [a'^A]A^ 

The  SVD  of  A  is  given  by 

A  =  ULV^ 


p.\GE3or:3 


(6) 


In  tenns  of  the  SVD  of  A  then,  (he  Moore-Penrose  inverse  is  given  by 

A*  -  VE+i/^ 

r/^nihmina  (1),  (4).  (5)  and  (6)  the  projection  operator  P  can  now  be  cast  in  terms  of  the 
singular  vectors  of  A  as 

Pm  AA*  UU^ 

The  sin<Tilar  vectors  which  make  up  the  columns  of  U  are  an  orthonoimal  basis  set  which 
span  the  range  of  A.  However,  not  all  of  the  singular  vectors  are  equaUy  important  m  the 
reconstruction  (projection)  process.  Those  associated  with  the  largest  singular  values  are 
the  most  important,  in  the  sense  that  excluding  them  would  increase  the  mean  square  error 
between  any  projected  data  vector  and  the  original  data  vector  Xj  by  the  l^est  amoimt. 
Conversely,  excluding  the  singular  values  associated  with  the  smallest  singular  values 
would  have  the  smallest  impact  on  the  mean  square  error.  Since  it  is  assumed  that  it  is  the 
interferers  which  contribute  the  most  to  the  structure  of  the  data,  the  singular  values  asso¬ 
ciated  with  their  presence  should  be  the  largest.  Therefore  we  construct  P  only  out  of  the 
singular  vectors  in  U  which  are  associated  with  the  lar^t  singular  values. 

Fi<mre  2  shows  a  power  spectral  density  plot  of  signal  0  constructed  from  2048  samples 
(approximately  1  second  of  data)  taken  form  a  starting  position  near  the  15  second  pomt  of 
fi«Tire  1.  The  presence  of  the  interfering  tones  is  clearly  visible,  and  it  is  noted  that  they  are 
not  all  harmonically  related.  Note  also  that  the  bulk  of  the  noise  floor  lies  below  about  -65 
dB  and  that  the  interferers  lie  primarily  in  that  region. 

Figure  3  shows  the  singular  values  of  a  data  matrix  A  constructed  from  a  set  of  121  sample 
v^tors  {x  }  ralfpn  from  signal  0.  Each  sample  vector  was  2048  sample  points  in  length, 
and  there  was  a  50  percent  data  overlap  between  adjacent  vectors.  We  can  see  that  the 
power  of  the  associated  singular  vectors  extends  down  to  about  where  we  would  like  to 

Power  Spectral  Density  Of  Unprocessed  Silence  (Signal  0) 


PAGn  4  OF  23 


go.  so  we  will  use  all  of  the  singular  values  (the  A  matrix  is  of  rank  121)  to  compute  the 
projector  P  from  (7). 

However,  it  is  not  the  projection  of  the  interferers  which  we  want,  it  is  their  elimination. 
To  do  this  we  construct  the  null  space  projector  P'^.  This  projector  projects  a  given  vector 
into  the  null  space  of  A.  In  the  case  of  a  raw  data  vector  Xj  extracted  from  signal  0,  the 
result  of  the  projection  -  P^Xj  will  be  the^noise  (sUence)  only,  since  the  interferers  are 
in  the  range  of  A.  The  null  space  projector  P^  is  given  by 


where  I  is  the  identity  matrix. 

Figure  4  shows  the  spectrogram  resulting  from  the  null  space  projection  of  the  data 
graphed  in  figure  1.  Note  that  the  absolute  difference  in  grey  scale  is  not  as  important  as 
the  relative  difference.  Comparison  of  figure  4  with  the  original  spectrogram  of  figure  1 
shows  that  the  strongest  interferers  have  been  suppressed  the  most  while  the  timing  marks 
and  the  short  segment  of  signal  8  (300  Hz  tone)  above  80  seconds  have  been  retained  smce 
they  are  in  the  null  space  of  A. 

Figure  5  shows  the  corresponding  power  spectral  density  plot.  The  power  in  figure  5  is  not 
normalized,  however  the  range  is  the  same  in  order  to  facilitate  direct  comparison  vdth  fig- 
ure.2.  Note  that  the  interfering  tones  are  gone,  except  for  those  which  were  approximately 
60  dB  down  from  the  peak  power  (0  dB)  in  figure  2.  Note  also  that  apart  from  the  suppres¬ 
sion  of  the  frequencies  near  DC.  the  overaU  shape  of  the  noise  fioor  has  been  preserved. 
The  suppression  of  the  low  frequency  terms  is  a  consequence  of  usmg  ail  of  the  largest 
121  singular  vectors  to  construct  the  projector  matrix,  rather  than  being  more  selective. 
Future  work  will  address  this  issue.  There  is.  for  example,  a  marked  knee  in  die  curve  of 
figure  3  after  the  first  singular  value.  This  implies  that  the  first  singular  value  is  related  to 
the  DC  and/or  lower  frequency  terms  in  figure  2.  Using  the  singular  vector  corresponding 
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Processed  Signal  0 


Figure  4:  Spectrogram  of  processed  silence. 


to  this  tenn  in  the  construction  of  the  projection  operator  P  in  future  studies  may  prevent 
the  suppression  of  the  lower  frequencies. 


Power  Spectral  Density  Of  Processed  Silence  (Signal  0) 
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Figiire  5:  Power  spectral  density  of  processed  signal  0. 
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Test  Results 


This  outlines  the  test  results  to  date.  The  null  space  projector  of  the  prewous  sec¬ 

tion  is  used  to  eliminate  the  interferers  from  flie  signal  list  shown  m  table  2.  Each  of  these 
cases  will  be  considered  separately.  In  each  case,  a  spectrogram  and  power  spectral  den¬ 
sity  plot  of  the  original  signal  will  be  shown.  foUowed  by  similar  plots  for  the  processed 

signal. 
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Signal  18  -  A  multi-tone  signal  consisting  of  100  and  500  Hz  tones. 

Signal  18 : 100  And  500  Hz  Tones 
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Figure  6:  Spectrogram  of  unprocessed  signal  18 


Figure  7:  Power  spectnil  density  of  unprocessed  signal  18 


Processed  Signal  18 : 100  And  500  Hz  Tones 
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Figure  8:  Spectrogram  of  processed  signal  18 


Power  Spectral  Density  Of  Processed  Signal  18  ;  100  And  500  Hz  Tones 
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Figure  9:  Power  spectral  density  of  processetl  signal  18 
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Figure  10:  Spectrogram  of  unprocessed  signal  19 
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Figure  11:  Power  spectral  density  ot  unprocessed  signal  19 
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Figure  12:  Spectrogram  of  processed  signal  19 
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Figiire  13:  Power  spectral  density  of  processed  signal  19 
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lal  20  -  An  impulse  train  with  a  30  Hz  repetition  rate. 

Signal  20  : 30  Hz  Impulse  Train 


1000 


800 


600 


400 


200 


”0  10  20  30  40  50  60  70  c 

Time  (Seconds) 

Figure  14:  Spectrogram  of  unprocessed  signal  20 


Power  Spectral  Density  Of  Signal  20  : 30  Hz  Impulse  Train 
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Figure  15:  Power  spectral  density  of  unpn)cessed  signal  20 
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Processed  Signal  20  :  30  Hz  Impulse  Train 


Figure  16:  Spectrogram  of  processed  signal  20 
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Figure  17:  Power  spectral  density  of  processed  signal  20 
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Processed  Signal  21  : 65  Hz  Impulse  Train 
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Figtire  22:  Power  S!)ecfral  density  of  processed  signal  21 
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Signal  22  -  An  multi-impulse  train  with  a  15  and  65  Hz  repetition  rate. 

Signal  22  : 15  and  65  Hz  Multiple  Impulse  Train 
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Figure  23 1  Spectrognini  of  unprocessed  signui 
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Power  Spectral  Density  Of  Signal  22  : 15  And  65  Hz  Impulse  Trains 
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Figure  24:  Power  spectral  density  ofunpnKessed  signal  22 
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Figure  25:  Spectrogram  of  processed  signal  22 


Power  Spectral  Density  Of  Processed  Signal  22  : 15  And  65  Hz  Impulse  Trains 
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Signal  23  -  An  multi-impulse  train  witli  a  30  and  65  Hz  repetition  rate. 


Figure  27:  Spectrogram  of  unprocessed  signal  23 
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Processed  Signal  23  : 30  and  65  Hz  Multiple  Impulse  Train 
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Power  spectral  Density  Of  Processed  Signal  23  : 30  And  65  Hz  Impulse  Trains 
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Figure  30:  Power  spectral  density  of  processed  signal  23 


Signal  24  -  A  QPSK  signal  with  a  baud,  rate  of  50  Hz  centered  on  250  Hz 
Sional  24 :  QPSK  With  50  Hz  Baud  Rate  Centered  On  250  Hz 
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Figure  31:  Spectrogram  of  unprocessed  signal  24 
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Figure  32:  Power  spectral  density  of  unprocessed  signal  24 
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4 :  QPSK  With  50  Hz  Baud  Rate  Centered  On  250  Hz 


Figure  33:  Spectrogmiii  of  processed  signal  24 
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Figure  34:  Power  spectral  density  of  processed  signal  24 
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Figure  35:  Spectrogram  of  unprocessed  signal  2? 


Figure  36:  Power  spectr;il  density  of  unprocessed  signal  25 


Processed  Signal  25  :  QPSK  With  100  Hz  Baud  Rate  Centered  On  250  Hz 
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SPCOT 

MEMO: 

September  1, 1993 

Edward  C,  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.Real 

885-5147 

RE: 

piiminarinn  of  interfering  signals  from  acoustic  data  using  LMS. 

FAX  885-0631 

CC: 

This  report  outlines  the  ability  of  the  LMS  algorithm  to  capture  sisals  taken  from  the 
ARO  data  set  The  basic  operational  parameters  for  this  work  are  outlined  in  the  following 

table.  _ _ _ _ _ _ 


Parameter 

Value 

Comments 

Sample  Rate 

2083.33  Hz 

As  given  in  the  ARO  test  plan. 

Data  Acquisi¬ 
tion  Date 

Dec.  7. 1992 

As  given  on  data  tape  label. 

Time 

1  to  1:10  pm 

ibid. 

Location 

North  Tower. 

WSMR 

ibid. 

Tape  #  That 
Data  Is  On 

2 

ibid. 

Qiannel  Used 

3 

Relatively  clean  data. 

Group  #  Used 

1 

There  are  2  groups  on  this  tape. 

Table  1:  Basic  operational  parameters 


For  this  study  the  signal  extracted  G>y  hand)  from  this  data  was  signal  0  (silence).  The 
LMS  algorithm  was  taken  from  “Adaptive  Filters”  by  CJ.N.  Cowan  and  P.M.  Grant. 
Prentice-Hall  Signal  Processing  Series.  Alan  V.  Oppenheim.  Series  Editor,  1985.  The 


Parameter 

Value 

Comments 

Number  Of 
Taps 

15 

As  directed. 

Learning 
Constant  (u) 

3.7509e-ll 

Initial  estimate.  Equal  to  ^ —  where 

^  max 

X  is  the  largest  eigenvalue  of  a  data 

matrix  composed  of  121  sample  vectors 

of  2048  points  each  with  50  percent 

overlap,  taken  from  a  segment  of  the 

signal.  This  value  was  adjusted  upwards 

as  is  shown  below.  See  above  reference 

for  rational  for  choosing  u  this  way. 

1 

fable  2:  Parameters  of  the  LMS  algorithm 

P.\GE  1  OF9 


Preliminaries  -  The  LMS  Algorithm 

The  update  scheme  for  the  N  element  coefficient  vector  h  obtmed  via  the  adaptive  LMS 
algoridim  is  given  in  the  above  reference  and  is  repeated  here  for  convemence. 

h(n)  h(n-l)  +2ue^ 

fc(n-l) 

where  n  is  the  iteration  number,  u  is  the  learning  constant  mentioned  above,  e,  is  the 
error  term  for  the  n*  iteration  given  as 

(n)  •=  yin) (n)h(n) 

A(«) 

and  where  v(n)  is  the  data  point  and  $  («)  is  an  N  element  data  vector  consisting  of 


and  where  y  (n)  is  the  n 
the  following  data  points 


y(n-l) 
yin -2) 

yin-N)_ 


Results 

Employing  this  algorithm,  and  training  on  10  seconds  of  silence  with  a  lea^g  constant  u 
of  3.75W-II  resulted  in  the  error  curve  shown  m  figure  1.  Figure  2  shows  the  corre- 

spending  coefficient  vector  h . 


Error  curve  for  a  learning  constant  of  3.7509e-1 1 


3 

S  600 

(D 


Iteration 


Fi«ure  1:  Error  curve  for  a  learning  constant  of  3.7509e-ll 
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Filter  coeffidents  for  a  learning  constant  of  3.7509o-1 1 


Figure  2:  Coefficient  vector  for  a  leamiiig  constant  of  3.7509e-ll. 


Fionres  3  and  4  show  similar  graphs  for  a  learning  constant  of  1.2003e-9  (32  *  u)  wble 
fi^es  5  and  6  show  graphs  for  a  learning  constant  of  4.8011e-9  (128  *  u).  For  clarity,  fig¬ 
ure  7  shows  a  blow  up  of  the  error  curve  plot  of  figure  5. 

Figures  8  and  9  shows  graphs  for  a  learning  constant  of  9.6022e-09  (256  *  u),  note  the 
change  in  scale  on  figure  8. 

Error  curve  for  a  learning  constant  of  1 .2003e-09 
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Figure  3:  Error  curve  for  a  leaniing  constant  of  12003e-9 
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Fflter  coeHidents  for  a  learning  constant  of  1  -20036-9 
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Figure  4:  Coefficient  vector  for  a  learning  constant  of  12003e-9 

Figures  10  and  11  show  the  graphs  for  a  learning  constant  of  1.9204e-08  (512  *  u). 

Figures  12  and  13  are  for  a  learning  constant  of  3.8409e-08  (1024  *  u).  Note  that  at  this 
learning  constant  is  driving  the  algorithm  into  instability.  Figure  14  shows  a  blow 

up  of  figure  12. 


Error  curve  for  a  learning  constant  of  4.8011  e-9 
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Fi«ure  5:  Error  curve  for  a  learning  constant  of  4.80lle-9 
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FBter  coefficients  for  a  learning  constant  of  4.801 1o-9 


Error  Value  (e(n)) 


Filter  coefficients  for  a  learning  consfant  of  9.6022e-09 


Figure  9:  Coefficient  vector  for  a  learning  constant  of  9.6022e-9 
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Error  Value  (e(n)) 


Error  curve  for  a  learning  constant  of  1 .9204e*08 


Filter  coefficients  for  a  learning  constant  of  1.9204e-08 
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MEMO:  September  17. 1993 


Edward  C.  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.Real 

885-5147 

\  RE: 

Plots  from  file  lO.OS.OO.gl  of  the  NASA  wind  tunnel  data. 

FAX  885-0631 

CC: 

S.Lang 

1 

I 

I 

I 

I 


Hiis  memo  contniiis  the  icQuested  plots  from  the  NASA  wind  tunnel  data  outlined  in  the 
September  15  fax  to  Doug  Mook  from  Henry  Abarbanel  of  UC  San  Diego.  The  specifics 
of  the  data  plotted  are  contained  in  the  table  below. _ 


Parameter 

Value 

Comments 

Sample  Rate 

2083.33  Ifr 

As  given  in  the  fax  mentioned 
above. 

File  Name 

10.05.00.gl 

From  NASA  wind  tunnel  data 
accompanying  above  fax. 

Data  Creation 
Date 

July  9. 1993 

As  specified  in  the  extracted  tar 
file  date  stamp. 

Data  Creation 
Time 

15:26 

ibid. 

Location  Of  Sen¬ 
sor  Array 

Front  center 
of  canopy 

6"  diameter  circular  sensor 
array  (see  fax). 

Sensor  Number 
Within  Array 

4 

As  specified  by  C.  Myers.  There 
are  8  sensors  total. 

Channel  Used 
For  Sensor  Output 

ibid.  There  are  8  channels  total. 

Group  #  Used 

1 

ibid. 

Table  1:  Data  Specifics 


Figure  1  is  a  spectrogram  of  the  entire  data  stream  outlined  in  table  1.  Figure  2  is  a  power 
spectral  density  plot  of  the  same  data  over  its  entire  length  (125952  point  DFT). 


Ed. 
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Frequency  (Hz) 


File  10.05.Ql  :  Sensor  4  (Channel  3) 


Time  (Seconds) 

Figure  1:  Spectrogram  of  channel  3  of  data  file  10.05.00.gl 


File  10.05.gl  ;  Sensor  4  (Channel  3) 
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Figure  2:  PSD  plot  of  channel  3  of  data  file  10.05.00.gl 
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SPCOT 

MEMO: 

September  19, 1993 

Edward  C.  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.Real 

Plots  from  file  10.05.00.gl  of  the  NASA  wind  tunnel  data  -  all  channels. 

885-5147 

RE: 

FAX  885-0631 

CC: 

S.Lang 

Tlis  „«mo  con,3i»,  plo«  f-o»  cba»n.U  of  wM 

oiiaincd  in  Ihe  Sepranber  15  fai  to  Doug  Mook  ftom  Heniy  Abatbanel  of  UC  San  Diego. 


1  Parameter  1 

Value 

Comments 

2083.33  Hz 

As  given  in  the  fax  mentioned 
above. 

File  Name 

10.05.00.gl 

From  NASA  wind  tunnel  data 
accompanying  above  fax. 

Data  Creation 
Date 

July  9, 1993 

As  specified  in  the  exuacted  tar 
file  date  stamp. 

15:26 

ibid. 

Location  Of  Sen¬ 
sor  Array 

Front  center 
of  canopy 

6”  diameter  circular  sensor 
array  (see  fax). 

Channels  Used 
For  Sensor  Output 

aU 

As  specified  by  C.  Myers.  There 
are  8  channels  total. 

Group  #  Used 

1 

ibid. 

Ficure  1  3  5  7  9, 11, 13,  and  15  are  spectrograms  of  the  entire  data  streams  for  each  of 
"^'::lcis  ouaine6  in  table  1.  Figure  2. 4, 6,  8, 10, 12. 14.  and  16  ^e  cormsgn^ng 
power  spectral  density  plots  of  the  same  data  (over  the  enure  data  length  -  1.5952  point 

npro 


Ed. 
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Figure  8:  Power  spectral  density  plot  of  channel  4  of  file  lO.OS.OO.gl 
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Magnitude  (dB) 


RIe  10.05.00.g1  Channel  5 


Figure  9:  Spectrogram  of  channel  5  of  file  10.05.00.gl 


File  10.05.00.g1  Channel  5 


Figure  10:  Power  spectral  density  plot  of  channel  5  of  file  10.05.00. 
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SPCOT 

MEMO: 

October  1, 1993 

Edward  C.  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.Real 

Plots  from  file  10.05.00.gl  of  the  NASA  wind  tunnel  data  -  Oiannel  7 

885-5147 

RE: 

FAX  885-0631 

CC: 

S.Lang 

TOs  memo  contains  PSD  plots  from  channel  7  of  the  NASA  vrind  ^1  data^tlined  in 
the  September  15  fax  to  Doug  Mook  from  Henry  Abarbanel  of  UC  San  Diego.  The  specif- 


JL  UUW  uaua  - - 

Parameter  Value 

Comments 

Name 

10.05.00.gl 

From  NASA  wind  tunnel  data 
accompanying  above  fax. 

Data  Creation  Date 

July  9. 1993 

As  specified  in  the  extracted  tar  file 
date  stamp. 

Data  Creation  Time 

15:26 

ibid. 

Location  Of  Sensor 
Array 

Front  center 
of  canopy 

6*’  diameter  circular  sensor  array 
(see  fax). 

Channel  Used  For 
Sensor  Output 

7 

As  specified  by  C.  Myers.  There  are 

8  channels  total. 

Group  #  Used 

1 

ibid. 

Sample  Rate 

2083.33  Hz 

As  given  in  the  fax  mentioned 
above. 

EFT  size  used  in  PSD 
computation 

2048 

Approximately  1  Hz  resolution,  and 
approximately  1  second  of  data. 

Starting  Data  Point 

1 

Bejunning  of  data  set. 

FIR  filter  order  used 

15 

•  •  a  ^  A. 

All  prediction  error  plots. 

-  o _ 

Overview: 


A  15  tap  FIR  LS  predictor  filter  was  created  in  madab  (file  Isfilterjn)  and  used  to  process 
the  <’hannol  data  for  a  variety  of  vector  (d)  and  predictor  (k)  shifts.  These  values  are  out- 
lined  ia  table  2. 


Vector  Shift 
(d) 

Predictor  Shift 
(k) 

1 

1 

130 

130 

The  meaning  of  these  parameters  is  defined  in  the  Mowing  discussion. 


RACE  1  OF  7 


Where  x(n)  is  the  n-th  data  point,  d  is  the  vector  shift  parameter  (defines  the  time  shift 
between  column  vectors  in  X),  L  is  the  number  of  taps  in  the  FIR  ater  (15  m  the  cases 
presented  here),  a  is  the  tap  coefficient  vector,  and  k  is  the  predictor  shift  parameter 
(defines  how  far  into  the  future  that  the  prediction  is  to  be  nmde).  Tte  LS  coefficients  of 
the  filter  (a)  are  calculated  via  the  Moore-Penrose  inverse  as  in  equation  2. 

a  -  (X"X)"*X"x  (2) 

The  error  (e)  is  then  calculated  as  in  equation  3 

e  =  x-Xa  (3) 


RAGE  2  OF  7 


Case  1  -  Original  PSD; 

Hie  plot  presented  here  is  the  (unweighted)  power  spectral  density  (PSD)  of  the  first  2048 
points  of  the  data  from  channel  7.  This  was  done  to  supply  a  basis  for  companson  with  the 

other  plots. 


PSD  Of  File  10.05.00.g1  Channel  7 


Figure  1:  PSD  of  first  2048  data  points  of  channel  7  of  file  10.05.00.gl 
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Case  2  -  d  =  1,  k  =  1: 

For  this  case  a  vector  shift  (d)  and  a  predictor  shift  (k)  of  one  were  used.  Ilie  plot  here  and 
for  aU  subsequent  cases,  represents  the  PSD  of  the  first  2048  data  ^mts  of  the  resulting 
error  function  (equation  3).  The  length  of  the  column  vectors  of  X  (number  of  rows)  was 
60,000  samples  in  all  cases. 


PSD  Of  D  =  1  K  =  1  Error  For  File  10.05.00.g1  Channel  7 


Figure  2:  PSD  of  error  curve  for  parameters  d  -  1,  k  -  1. 
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Power  (dB) 


Case  4  -  d  =  130,  k  =  1 


PSD  0(  D  =  130  K  «=  1  Error  For  RIe  10.05.00.g1  Channel  7 


Figure  4:  PSD  of  error  curve  for  parameters  d  =  130,  k  =  1. 
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SPCOT 

MEMO: 

October  3, 1993 

Edward  C.  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.Real 

885-5147 

RE: 

Embedding  Himpnsinn  plots  of  residual  error  resulting  from  a  15  tap  LS  FIR 
filter. 

S.Lang 

FAX  885-0631 

_ a  PP* 

This  memo  contains  percentage  of  joint  false  neighbor  Vs.  embedding  dimension  plott. 
IMs  was  g«»np.rateH  from  residual  error  data  resulting  from  a  15  tap  least  square  (^) 
FIR  predictor  filter  (see  October  1, 1993  memo  to  C.  Myers  from  E.  Real).  The  original 
data  was  taken  from  channel  7  of  the  NASA  wind  tunnel  data  outlined  in  the  September 
15  fax  to  Doug  Mook  from  Henry  Abarbanel  of  UC  San  Diego.  Tire  specifics  of  the  data 


Parameter 

Value 

Comments 

File  Name 

10.05.00.gl 

From  NASA  wind  tuimel  data 
accompanying  above  fax. 

Data  Creatioa  Date 

July  9, 1993 

Data  Creation  Time 

15:26 

ibid. 

Location  Of  Sensor 
Array 

Front  center 
of  canopy 

6”  diameter  circular  sensor  array 
(see  fax). 

rhflnnel  Used  For 
Sensor  Output 

7 

As  specified  by  C.  Myers.  There  are 

8  channels  total. 

Group  #  Used 

1 

ibid. 

Sample  Rate 

2083.33  Hz 

As  given  in  the  fax  mentioned 
above. 

FFT  size  used  in  PSD 
computation 

2048 

Approximately  1  Hz  resoludon.  and 
approximately  1  second  of  data. 

Starting  Data  Point 

1 

Beginning  of  data  set. 

FIR  filter  order  used 

15 

All  prediction  error  plots. 

Table  1:  Data  Specifics 
filter  used  and  the  meaning  of  parameters,  please  see  the  October  1. 1993  memo. 
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Percentage  Of  Joint  False  Neighbors 


Embedding  Dimension  Plot  for  d  =  1.  k  =  Is 


Figure  1:  Pereentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  1  Predictor  Shift  (k)  =  1 


RAGE  2  OF  5 


Embedding  Dimensioii  Plot  for  d  =  1,  k  =  130 


Parameters  d  - 1 .  k  - 130 :  Data  Residual  Error,  15  Tap  LS  FIR  Filter 


Figure  2:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  1  Predictor  Shift  (k)  =  130 
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Percentage  Of  Joint  False  Neighbors 


Embedding  Dimension  Plot  for  d  =  130,  k  =  1 


Figure  3:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  130  Predictor  Shift  (k)  =  1 
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Embedding  Dimension  Plot  for  d  =  130,  k  =  130 


Figure  4:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  130  Predictor  Shift  (k)  =  130 
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SPCOT 

MEMO: 

October  6. 1993 

Edward  C.  Real 

TO: 

C.  Myers 

PTP2-A001 

FROM: 

E.  Real 

PSD  plots  of  residual  error  resulting  from  a  1  tap  LS  FIR  filter. 

885-5147 

RE: 

FAX  885-0631 

CC: 

S.Lang 

This  memo  shows  residual  error  PSD  plots  resuUing  from  a  1  tap  le^t^uare  (LS)  FIR 
predictor  filter  (see  October  1,  1993  memo  to  C.  Myers  from  E  Real),  m  on^l  data 
was  taken  from  channel  7  of  the  NASA  wind  tunnel  data  outlined  m  the  September  15  fax 
to  Doug  Mook  from  Henry  Abarbanel  of  UC  San  Diego.  Die  specifics  of  the  data  plot«^ 


Parameter 

Value 

Comments 

File  Name 

10.05  .OO.gl 

From  NASA  wind  tunnel  data 
accompanying  above  fax. 

Data  Creation  Date 

July  9. 1993 

As  specified  in  the  extracted  tar  file 
date  stamp. 

Data  Creation  Time 

15:26 

ibid. 

Location  Of  Sensor 
Array 

Front  center 
of  canopy 

6"  diameter  circular  sensor  array 
(see  fax). 

Channel  Used  For 
Sensor  Output 

7 

As  specified  by  C.  Myers.  There  are 

8  channels  total. 

Group  #  Used 

1 

ibid. 

Sample  Rate 

2083.33  Hz 

As  given  in  the  fax  mentioned 
above. 

FFT  size  used  in  PSD 
computation 

2048 

Approximately  1  Hz  resolution,  and 
approximately  1  second  of  data. 

Starting  Data  Point 

1 

Beginning  of  data  set. 

FIR  filter  order  used 

1 

All  prediction  error  plots. 

c'  *c. _ 

used  and  the  m^nnino  of  parameters,  please  see  the  October  1. 1993  memo. 
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Power  (dB) 


PSD  plot  of  original  (unprocessed)  data: 


PSD  01  File  10.05.00.g1  Channel  7 


Figure  1:  PSD  of  data  from  file  10.05.00.gl.  channel  7 
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Power  (dB) 


PSD  Plot  for  d  =  1,  k  =  1 


PSD  Of  D  =  1  K  =  1  Error  For  File  10.05.00.g1  Channel  7 


Frequency  (Hz) 

Figure  2:  PSD  of  residual  error  resulting  from  a  1  tap  FIR  filter 
Vector  Shift  (d)  =  I  Predictor  Shift  (k)  =  1 


PSD  Plot  for  d  =  1,  k  =  130 
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PSD  Plot  for  d  =  130,  k  =  1 


PSD  Of  D  =  130  K  =  1  Error  For  File  10.05.00.g1  Channel  7 


Figure  4:  PSD  of  residual  error  resulting  from  a  1  tap  FIR  filter 
Vector  Shift  (d)  =  130  Predictor  Shift  (k)  =  1 


PSD  Plot  for  d  =  130.  k  =  130 


Figure  5:  PSD  of  residual  error  resulting  from  a  1  tap  FIR  filter 
Vector  Shift  (d)  =  130  Predictor  Shift  (k)  =  130 
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SPCOT 


MEMO: 


October  7, 1993 


Edward  C.  Real 


TO:  C.  Myers 


PTP2-A001 


885-5147 


FAX  885-0631 


FROM:  E.  Real 

RE:  Embedding  dimension  plots  of  residual  error  resulting  from  a  1  tap  LS  FIR 

filter. 

CC:  S.Lang.W.Nfistretta  _ 


This  memo  contains  percentage  of  joint  false  neighbor  Vs.  embedding  dimension  plo« 
This  data  was  generated  from  residual  error  data  resulting  from  a  1  tap  least^uare  (LS) 
FIR  predictor  filter  (see  October  1.  1993  memo  to  C.  Myers  from  E.  Real).  Tfre  ongi^ 
data  was  taken  from  channel  7  of  the  NASA  wind  tunnel  data  outhned  m  die  Sepffiinber 
15  fax  to  Doug  Mook  from  Henry  Abarbanel  of  UC  San  Diego.  The  specifics  of  the  data 


Parameter 

VaiiiP  Comments 

File  Name 

10.05.00.gl 

From  NASA  wind  tunnel  data 
accompanyum  above  fax. 

Data  Creadon  Date 

July  9. 1993 

As  specified  in  the  extrxted  tar  file 
date  stamp. 

Data  Creadon  Tune 

15:26 

ibid. 

Locadon  Of  Sensor 
Arrav 

Front  center 
of  canopy 

6"  diameter  circular  sensor  array 
(see  fax). 

Channel  Used  For  | 
Sensor  Output 

7 

As  specified  by  C.  Myers.  There  are 

8  channels  total. 

Group  #  Used 

1 

ibid. 

Sample  Rate 

2083.35  Hz 

As  given  in  the  fax  mendoned 
above. 

FFT  size  used  in  PSD 
computadon 

2048 

.Approximately  1  Hz  resoludon.  and 
approximately  1  second  of  data. 

Starting  Data  Point 

1 

Beginning  of  data  set. 

FIR  filter  order  used 

1 

.All  piedicdon  error  plots. 

r*  _  _  ’C _ 

filter  used  and  the  meaning  of  parameters,  please  see  the  October  1, 1993  memo. 

Note  that  for  the  single  tap  filter  the  vector  shift  parameter  (d)  is  irrelevant  since  the  data 
matrix  contains  only  a  single  column.  However,  plots  for  the  d  =  1  and  d  =  130  case  are 
inciuded  for  compieteness. 


i>aohioi*.s 


Percentage  Of  Joint  False  Neighbors 


Embedding  Dimension  Plot  for  d  =  1,  k  =  1: 


Parameters  d  -  1.  k  =  1  :  Data  Residual  Error.  1  Tap  LS  FIR  RIter 


lOh . 

! 

°0  2  4  6  8  io  12  U  Te 

Embedding  Dimension 

Figure  1:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  1  Predictor  Shift  (k)  =  1 
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Percentage  Of  Joint  False  Neighbors 


Fmh<»HHing  Dimension  Plot  for  d  =  1,  k  =  130 


Parametere  d  =  1,  k  =  130 :  Data  Residual  Eiror.  1  Tap  LS  FIR  Filter 
100| - ^ - ' - ' - ! - -  ^  ’ 


loi- 


2  4  '  i  ^  8  ib  l'2  ^4.  1^6 

Embedding  Dimension 

Figure  2:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  I  Predictor  Shift  (k)  =  130 
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Embedding  Dimension  Plot  for  d  =  130,  k  -  1 
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Figure  3:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  130  Predictor  Shift  (k)  =  1 
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Embedding  Dimension  Plot  for  d  =  130,  k  =  130 


Parameters  d  =  130,  k  =  130 :  Data  Residual  Error,  1  Tap  LS  FIR  Fille  r 
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Figure  4:  Percentage  Of  Joint  False  Neighbors  Vs.  Embedding  Dimension 
Vector  Shift  (d)  =  130  Predictor  Shift  (k)  =  130 
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